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1 Introduction 


Orthonormal polynomials on the real line are defined by the orthogonality conditions 


Pn(x)p m (x) dp{x) = 5 m ,n, 771, 71 > 0, 


( 1 . 1 ) 


where /i is a positive measure on the real line for which all the moments exist and p n (x) = 
r y n x n + • • •, with positive leading coefficient > 0. A family of orthonormal polynomials 
always satisfies a three-term recurrence relation of the form 


xp n {x ) = a n+ ip n+ i(x) + b n p n (x ) + a n p n -i(x), n > 0 , 

with p_i = 0 and 

Po = 7o = ^ J dfj,(x) 

Comparing the leading coefficients in (1.2) gives 


- 1/2 


®n+1 


In 


I'n+l 


> 0 , 


and computing the Fourier coefficients of xp n (x) in (1.2) gives 


dri. 


b n — 


XPn(x)Pn-l(x) dh(x), 

xpl(x)dp(x). 


( 1 . 2 ) 


(1.3) 


(1.4) 


(1.5) 

( 1 . 6 ) 


The converse statement is also true and is known as the spectral theorem for orthogonal 
polynomials', if a family of polynomials satisfies a three-term recurrence relation of the 
form (1.2), with a n > 0 and b n E M and initial conditions p 0 = 1 and p_i = 0, then there 
exists a probability measure /i on the real line such that these polynomials are orthonormal 
polynomials satisfying (1.1). This gives rise to two important problems: 
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Problem 1. Suppose the measure // is known. What can be said about the recurrence 
coefficients {a n : n — 1, 2, 3,...} and {b n : n — 0,1, 2,...}? This is known as the 
direct problem, for orthogonal polynomials. 


Problem 2. Suppose the recurrence coefficients {a n+ i,b n : n = 0,1,2,...} are known. 
What can be said about the orthogonality measure pi This is known as the inverse 
problem, for orthogonal polynomials. 


The recurrence coefficients are usually collected in a tridiagonal matrix of the form 


/b 0 a i \ 

o. i b\ ct 2 

Oj2 &2 «3 

I a 3 b 3 

\ ■■■/ 


(1.7) 


which acts as an operator on (a subset of) £2 (N) and which is known as the Jacobi matrix 
or Jacobi operator. If J is self-adjoint, then the spectral measure for J is precisely the 
orthogonality measure /j. Hence problem 1 corresponds to the inverse problem for the 
Jacobi matrix J and problem 2 corresponds to the direct problem for J. 

In the present paper we will study problem 1 for a few special cases. In Section 
2 we study measures on the real line with an exponential weight function of the form 
dp{x) = \x\ p exp(— \x\ m ) dx, which are known as Freud weights, named after Geza Freud 
who studied them in the 1970’s. It will be shown that the recurrence coefficients a n 
satisfy a non-linear recurrence relation which corresponds to the discrete Painleve I equa¬ 
tion and its hierarchy. In Section 3 we will study a family of orthogonal polynomials on 
the unit circle. We will first give some background on orthogonal polynomials on the 
unit circle and the corresponding recurrence relations. We will study the weight func¬ 
tion w{6) = exp(Acosd), and it will be shown that the recurrence coefficients satisfy a 
non-linear recurrence relation which corresponds to discrete Painleve II. These orthog¬ 
onal polynomials play an important role in random unitary matrices and combinatorial 
problems for random permutations. In Section 4 we will study certain discrete orthogonal 
polynomials related to Charlier polynomials. The recurrence coefficients a n and b n are 
shown to satisfy a system of non-linear recurrence relations which are again related to 
the discrete Painleve II equation. Finally, in Section 5 we consider certain (/-orthogonal 
polynomials which are (/-analogs of the Freud weight. We will show that the recurrence 
coefficients satisfy a g-deformed Painleve I equation. Most of the material in this paper is 
not new: the recurrence relations in Section 2 were already obtained by Freud in [7] and 
are known as Freud equations in the held of orthogonal polynomials. It was Magnus [17] 
who made the connection with the discrete Painleve equation I. The recurrence relation in 
Section 3 was found by Periwal and Shevitz [20] (see also Hisakado [9], Tracy and Widom 
[23]; Baik [1] used the Riemann-Hilbert approach to obtain the Painleve equation). The 
recurrence relations in Section 4 were obtained by Van Assche and Foupouagnigni in 
[25]. The results in Section 5 were first obtained by Nijhoff [19]. We hope that bringing 
together these various orthogonal polynomials and the corresponding discrete Painleve 
equations will be illuminating and encourage researchers in the held of orthogonal poly¬ 
nomials and researchers in integrable difference equations to talk to each other, and that 
the interaction between both areas of mathematics will shed some extra light on either 
subject. 
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2 Freud weights 

Freud weights are exponential weights on the real line (—oo, oo) of the form 

w p (x) = |a:| p exp(— |a:| ,ri ), p > —1 pm > 0. 

They were considered by Geza Freud in his 1976 paper [7], where he gave a recurrence 
relation for the recurrence coefficients a n when m = 2,4,6. For these cases Freud found 
the asymptotic behavior of the recurrence coefficients a n and he formulated a conjecture 
for this asymptotic behavior for every m > 0. This conjecture really started the analysis 
of general orthogonal polynomials on unbounded sets of the real line, since before Freud’s 
work only very special cases such as the Hermite polynomials and the Laguerre polynomi¬ 
als were studied in detail. In this section I will repeat Freud’s analysis of the recurrence 
coefficients of Freud weights, make the connection with discrete Painleve equations (which 
was not known to Freud but first pointed out by Magnus in [17]), and point out what has 
been done after Freud. Observe that Freud weights are symmetric, i.e., w p (—x) = w p (x), 
which implies that b n = 0 for n > 0. 


2.1 Generalized Hermite polynomials 


The case rn — 2 corresponds to generalized Hermite polynomials (and p = 0 are the 
Hermite polynomials). Generalized Hermite polynomials were already investigated by 
Chihara in [3] (see also [4]). The weight w 0 (x) = exp(— x 2 ) satisfies the first order differ¬ 
ential equation 

[w 0 (x)]' = -2 xw 0 (x) (2.1) 

which is the Pearson equation for the Hermite weight. In general a weight w satisfying 
a Pearson equation of the form \a{x)w{x)\' = r(x)w(x), with a a polynomial of degree 
at most two and r a polynomial of degree one, is called a classical weight. The weight 
functions for Hermite polynomials, Laguerre polynomials, and Jacobi polynomials are 
the classical weights for o of degree zero, one and two, respectively. Bessel polynomials 
appear for cr(x) = x 2 , but they are not orthogonal on the real line with respect to a 
positive measure. Freud’s idea was to compute the integral 


(Pn(x)p„-l(x)\x\ 1 ')' W 0 {x)dx 


( 2 . 2 ) 


in two different ways. The first way is simply working out the derivative in the integrand 
and to use the orthogonality to evaluate the resulting terms. This gives 

POO POO 


(Pn(UP»-lMM P ) , '“’oU) dx = / p , n (x)p n - 1 (x)\x\ P W 0 (x)dx 


(2,3) 


+ J 

+ p 


—oo 
oo 


Pn(x)p' n _ 1 (x)\x\ P Wo(x) dx 
p n {x)p n - i(x) 


X 


I^IGuo^) dx. 


(2.4) 

(2.5) 


For the right hand side in (2.3) we use the fact that 


Pn( X ) 


n^ n x n 1 + lower order terms 
n ^- p n _ i(x) + lower degree terms. 

7n—1 
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This gives 


p' n (x)p n -i(x)\x\ p w 0 (x ) dx = n 


In 


Tn—1 


For the integral in (2.4) we see that p' n _ 1 (x) is a polynomial of degree n — 2 and hence by 
orthogonality the integral vanishes. For the integral in (2.5) we use the fact that the weight 
function w p is even, i.e., w p (—x) = w p (x), which implies that p n {—x) = (—1 ) n p n (x). This 
means that p n (x)/x is a polynomial of degree n — 1 when n is odd and p n _ \{x)/x is a 
polynomial of degree n — 2 when n is even. Hence when n is even the integral in (2.5) 
vanishes, and when n is odd we have 


Pn(x) 

X 


p n -\(x) + lower degree terms, 

7n-l 


so that 


where 


P r = p j^ An 


X 


A„ = 


7n-l 


0 if n is even, 


1 if n is odd. 

Combining these results and using the expression a n = 7 n -i/7n gives 

/ (p„(x)p„- 1 (x)\x\ 1 ’)'w 0 (x) dx = - + . 


( 2 . 6 ) 


Observe that this holds whenever w is a symmetric weight on the real line. A second way 
to compute the integral in (2.2) is to use integration by parts, combined with Pearson’s 
equation (2.1) for the weight. This gives 


(Pn(x)pn-l(x)\x\ p )'W 0 (x) dx = - Pn{x)p n -l(x)\x\ p Vj' 0 (x) dx 


—oo 

POO 


2 / xp n (x)p n _ 1 (x)\x\ p w 0 (x) dx 

J —OO 


(2.7) 


Combining (2.6) and (2.7) then gives 


2 n + pA r 
° n ~ 2 


( 2 . 8 ) 


so that a n = \Jn + pA n / a/ 2. For p = 0, which are the Hcrmite polynomials, one has 
a n = \Jn/2. For generalized Hermite polynomials one has 


lim 

oo 




\fn 


1 

7!' 


(2.9) 
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2.2 Freud weight m = 4 

Let us now consider the weight w p (x) = |a;| p exp(—x 4 ), where p > —1. The Pearson 
equation for the weight w o(x) = exp(—a; 4 ) is 

M*)]' = ~4x 3 w 0 (x), 


which is a first order linear differential equation with polynomial coefficients, but since 
the polynomial coefficient is now of degree 3, the weight Wq is no longer a classical weight 
but a semi-classical one. The equation (2.6) remains valid for this Freud weight, but 
integration by parts give a different result. Indeed 

/ OO P OO 

(p n (x)p n -i{x)\x\ p )'w 0 (x) dx = - p n (x)p n - 1 (x)\x\ p w' 0 (x) dx 

-OO j —OO 

/ OO 

X 3 Pn{x)p n -i(x)\x\ P Wo(x) dx. (2.10) 

-OO 

This integral can be computed by applying the three term recurrence (1.2), with b n = 0, 
repeatedly. Indeed 

x 2 p n (x) = a n+ 1 xp n+ 1 (x) + CbnXPn-rix) 

&n+l (ln+2Pn+2 ((£") T (u n _|_i T ^n)Vn{x') T Cl n (l n —\p n — 2 [x ), 


and 


X 3 p n {x) = a n+ 1 a n+ 2 xp n+ 2 (x) + (a 2 n+1 + a 2 n )xp n (x) + a n a n _ixp n _ 2 (x) 

®n+l®n+2®n+3Pn+3 (x) Qn+l(® n+2 ^n+1 + a n)Pn+l(^) 

T Un(® )l -|-l T T d n — l)Pn— 1 (x) T d n Pn— l®n— 2Pn— 3(*^)• 


From this one easily finds 


^ 3 Pn(^)Pn-i(^)kl P Wo(^) dx = a n (a 


n+l 


+ a n + a n-lJ 


( 2 . 11 ) 


This result holds for n > 1 if we define a o = 0. Observe that one can obtain this also by 
using the calculus of the Jacobi operator, since 


(■J 3 


X 3 Pn 


(x)p r 


\X)W n {X 


dx, 


and the quantity of interest is (J 3 ) n ,n- 1 - This computation is quite simple since it amounts 
to some simple matrix multiplications. Combining (2.6) and (2.10)-(2.11) then gives 

4 a n( a n + 1 + a n + a n- 1) = n + pA n - (2-12) 


This time we do not get a n explicitly, but instead we get a second order non-linear 
recurrence relation for the recurrence coefficients. The initial conditions are ao = 0 and 
for ai we require that pi(x) = xp 0 (x)/a 1 has norm one, which means 


7o 


— / x 2 w p {x) dx = 1, 
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and together with (1.3) this means 


Q>\ — 



1/2 


If we put x n = 2 a? n then clearly x 0 — 0, x n > 0 for n > 0 and 


■I'nX n + X n —f) Tl + pA n . 


(2.13) 


This recurrence relation is the discrete Painleve equation d-Pi 


T X n T X n —\ 


z n + 7 (-l) r 


Xn 


+ 


z n = an + /3, 


with a — 1, (3 — p/ 2 , 7 = —p /2 and 5 = 0, since we can write A n = (1 — (—l) n )/ 2 . 

The observation that Freud’s equation (2.12) is a discrete Painleve equation was not 
known to Freud but was pointed out much later by Magnus in [17]. This means that the 
equation has the discrete Painleve property, which is known as singularity confinement 
(see, e.g., [ 8 ]): 


Definition 2.1 (discrete Painleve property). If x n is such that it results in a sin¬ 
gularity for x n+ i, then there exists a p G N such that this singularity is confined to 
x n+ i ,..., x n+p . Furthermore x n+p+ i depends only on x n _i, x n _ 2 ,.... 


The usual Painleve property for differential equations is that the only movable sin¬ 
gularities (singularities which depend on the initial conditions) of solutions of a Painleve 
equation are poles. Poles are isolated singularities, hence a discrete version of poles as 
singularities is to require that singularities of a discrete equation are confined. This is the 
case for discrete Painleve I. Consider for instance d-Pi in the form 


*£n(*£n+1 T *£n T %n— l) 


then 


n 

Xn+1 X n X n —lj 

x n 


and if x n = 0 then we have a singularity for x n+ i which becomes ±oo. For x n+ 2 we then 
find = 1=00 and for x n+ 3 we have the indeterminate form (± 00 ) + (=|=oo). A more careful 
analysis is to put x n = e and to expand x n+ k in a Laurent series in e. This gives 


Xn 
^n+1 
Xn +2 

^n+3 
Xn+4 


e 

n 


Xn— 1 ^ 

n 71 +I ... o, 

-1~ %n— 1 H-C + C^(+) 

e n 

n + 3 


n 


e + 0(e 2 ) 


n 


n + 3 


x n -\ + 0(e). 


So we see that as e —» 0 the indeterminate form for x n+3 becomes 0, but it does not give 
a new singularity for x n+ 4 . The singularity is confined to x n+ i, x n+2 ,x n+3 . 
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The solution of d-Pi can not be obtained in a closed form, but one can say a few 
things about the behavior of the solution. Freud obtained the asymptotic behavior of the 
solution of (2.13) in the following way. Since x n > 0 for n > 0 we have 

x 2 n < x n (x n+1 + x n + x n -i) = n + pA„, 

so that x n j\/n is a bounded and positive sequence. Define A and B as the smallest and 
largest accumulation points 

0 < A = lim inf —A < lim sup —2= = B < oo. 

Tl^OO yjn ji —>oO 

Choose a subsequence n' such that x n '/\fn! —■► A, then as n' —> oo in (2.13) we have 


1 < A(2B + A). 


In a similar way we can choose a subsequence n" such that x n "/\/rf —> B, and as n' —> oo 
in (2.13) we then find 

B(2A + B) < 1. 

Together this gives 

B(2A + B) < A(2B + A), 

so that B 2 < A 2 . But since we already know that A < B, this implies that A = B and 
hence 

lim ^-2= = A = B. 

n—>oc y/n 

If we take the limit in the recurrence relation (2.13) then one finds 3 A 2 = 1 so that 
A = l/y/3. Recall that x n = 2 a 2 z , hence as a result we have 


lim 


n—xx) 77, 


1/4 


^ 12 ' 


(2.14) 


The recurrence relation (2.13), with initial conditions 


x 0 = 0, xi 


2F(a±e) 

r(i?) 


is very unstable for computing the recurrence coefficients. Lew and Quarles [13] showed 
that there is a unique positive solution of the recurrence relation (2.13) with Xo = 0. 
Hence a small error in x\ eventually destroys the positivity of x n . In Figure 1 we plotted 
the values of x n obtained from the recurrence relation (with p = 0) by using 30 digits 
accuracy. The x n are following the asymptotic behavior \Jn/3 quite well until n ~ 50 
and then large deviations from the true solution appear. 

Lew and Quarles proved the unicity by showing that there is an operator F acting 
on a Banach space of infinite sequences with xq = 0, such that the positive solution 
x = (xo, Xi, X 2 ,...) of (2.13) is a fixed point: F(x) = x. The operator F is then shown to 
be a contraction, so that the fixed point is unique. 

Observe that if we take a weight function of the form w p (x) = |a:| p exp(— x A + Xx 2 ), 
then the Pearson equation becomes 


[■iuo^)]' = (— Ax 3 + 2A x)w 0 (x), 
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Figure 1: The result of computing x n from d-Pi using 30 significant digits 


and a slight modification of the previous computations gives the recurrence 

4a 2 (a 2 +1 + a; + a 2 ^) — 2 A a; = n + pA n . 

If we put x n = 2a 2 y then the x n satisfy the discrete Painleve equation d-Pi with a = 1, 
f3 = p/ 2 , 7 = —p /2 and 5 = A. 


2.3 Freud weight m = 6 

For the Freud weight w p (x) = |a;| /: 'exp(—x 6 ) one can proceed in a very similar way. The 
Pearson equation now becomes 

[w 0 (x)]' = -6x 5 w 0 (x), 


so that 

/ OO POO 

(p n (x)pn-i(x)\x\ p y w 0 (x) dx = 6 / x 5 p n (a;)p n _i(x)|z| p wo(:r) dx. 

-oo J —OO 

The integral on the right is 


x 5 p n (a;)p n _i(a;)|a;| p u;o(a;) dx = (J 5 ) n , n _1 


and this is 


(J 5 


/n,n 


-i = a r 


/ 2 2 I 4 I q 2 2 

l°n-2 a n-l + °n-l + ^ a n-l a n 


, 2 2 
+ a n-l a n+l 
4 


a’ + 2a;a^ +1 + 


J 4 

4 n+l 


7 2 

hi+l L 


2 1 

n+2/ 




Together with (2.6) this gives 

/? 2 / 2 2 I 4 | rx 2 2 I 2 2 

Da nl a n-2 a rc-l+ °n-l ' ^ a n-l a n ' a n-l a n+l 

+ a* + 2a;a; +1 + af m + a 2 n+l a~ n+2 ) = n + pA„. (2.15) 

This is a fourth order non-linear recurrence relation for the recurrence coefficients. It is 
within the hierarchy of discrete Painleve I [6]. Freud was also able to obtain the asymptotic 
behavior for the a n in this case. Obviously 

6 a® < 6a^(a^_ 2 a^_ 1 +a^_ 1 +2a^_ 1 a^+a^_ 1 a^ +1 +a^+2a^a^ +1 +a^ +1 +a^ +1 a^ +2 ) = n+pA n , 
hence a n /n 1//6 is a positive and bounded sequence. If we define 

a 2 a 2 

0 < A = lim inf < lim sup = B < oo, 

/A 


1/3 


n-> OO n'/'J n‘ 

then by taking a subsequence that converges to A we find 

6A(5B 2 + 4 AB + A) > 1 
and by taking a subsequence that converges to B we ford 

6B(5A 2 + AAB + B 2 ) < 1. 

Combining both inequalities gives 

6B(5A 2 + 4 AB + B 2 ) < 6A(5B 2 + AAB + A) 
so that A 2 B + B 3 < AB 2 + A 3 . This is equivalent with 

(B - A)(A 2 +AB + B 2 ) <(B- A)AB. 

If we assume that A < B, then this would imply that A 2 + B 2 < 0, or A = B = 0, which 
is impossible (since A < B). Hence our assumption is false and A = B. Taking the limit 
in (2.15) gives 60H 3 = 1, hence 


lim 


n—>oo n 1 / 6 >^60' 


(2.16) 


2.4 Freud’s conjecture 


On the basis of (2.9), (2.14) and (2.16), Freud made the conjecture that for every m > 0 
and p > — 1 one has 


lim —j— 

n—xx) 77 ,V m 


/ r(m + 1 ) y 1/m 
v r (f)r(f + \)j 


(2.17) 


Furthermore Freud showed that if the limit exists for even m, then it is equal to the 
expression in (2.17). Freud could not prove the existence of the limit for even m > 8 
because the central coefficient a n in the non-linear recurrence relation does not occur 
sufficiently often and more non central coefficients a n ±k with k ^ 0 appear, making the 
recurrence no longer ‘diagonally dominant’. The simple trick using lim sup and lim inf 
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then no longer suffices to show that the limit exists. The proof of Freud’s conjecture for 
every even integer m was given by Alphonse Magnus [15] [16]. His proof still consists of 
obtaining a non-linear recurrence relation for the a n (the Freud equation, which is within 
the hierarchy of discrete Painleve 1), but a more subtle argument is used to prove the 
existence of the limit. Freud’s conjecture for general m > 0 was finally proved by Lubinsky, 
Mhaskar and Saff [14], The proof for general m > 0 no longer uses a recurrence relation 
for the recurrence coefficients but relies on the Mhaskar-Rakhmanov-Saff number and 
results of weighted polynomial approximation and an equilibrium problem of logarithmic 
potential theory with external held. 

For m an even positive integer, Mate, Neva! and Zaslavsky [18] obtained an asymptotic 
expansion of the form 

o oo 

a n _ ST' _£fc_ 

/ TYl ' j^2k ^ 

k=0 

where Co is the constant in Freud’s conjecture (2.17), but the other coefficients Ck with 
k > 0 are not explicitly known. Their analysis is again based on the non-linear recurrence 
relation for the recurrence coefficients. 


3 Orthogonal polynomials on the unit circle 

In this section we will consider orthogonal polynomials on the unit circle T = {z E C : 
|^| = 1}. A very good source for the general theory is the recent set of books by Barry 
Simon [22], The sequence of polynomials {<p n ,n = 0,1, 2,...} is orthonormal on the unit 
circle with respect to a weight w if 

2t r _ 

< Vn{z)ip m {z)w{0 ) dB = (3.1) 

These polynomials are unique if we agree to make the leading coefficient positive: 

^Pn (^) ^nZ T ' * * ? K'n 0. 

The monic polynomials are usually denoted by <F n (^) = tp n (z)/hi n . An important property, 
which replaces the three term recurrence relation for orthogonal polynomials on the real 
line, is the Szego recurrence 



Z$n(z) = ®n+l(z) + a n $* n (z), 


(3.2) 


where (z) = z n Q n (l/z) is the reversed polynomial and <F n is the polynomial <F n but with 
complex conjugated coefficients. In [22] the recurrence coefficients a n (n = 0,1, 2,3,...) 
are called Verblunsky coefficients. They are given by a n = — < h r)+ i(0) and they satisfy 
|ce n | < 1 for n > 0 and «_i = —1. An important relation between n n and a n was found 
by Szego: 

n 

K 2 n = X^ fe (°)| 2 ’ 

k=0 


from which it follows that 


K n +1 K n — | T^n+1 (0) | — ^n+ll^nl j 


10 



so that 




(3.3) 


3.1 Modified Bessel polynomials 

We will take a closer look at the orthogonal polynomials on the unit circle for the weight 

w(9) = exp(A cos 9), 0 G [— ir, n\. 

observe that w(—9) = w(9), which implies that the Verblunsky coefficients are real. 
Ismail [11, pp. 236-239] calls the resulting orthogonal polynomials the modified Bessel 
polynomials since the trigonometric moments of this weight are in terms of the modified 
Bessel function 

i /*2i r i /»7r 

— / e me w(9)d6=— / cos n# exp (A cos#) dO = I n (\). 

2vr Jo Jo 

These polynomials appear in the analysis of unitary random matrices [20, 23, 9] and play 
an important role in the asymptotic distribution of the length of the longest increasing 
subsequence of random permutations [2], 

Periwal and Shevitz [20] found a non-linear recurrence relation for the Verblunsky 
coefficients of these orthogonal polynomials (see also [9, 23]). If z — e ie then 


w(9) = exp f^(z + -)] := w(z), 


and this weight satisfies a Pearson equation of the form 


Consider the integral 


w\z) = ^(1 - ^)w(z). 


77- f Vn{z)Vn+\{z)w'(z)—, 


then by means of Pearson’s equation (3.4) we find 


(3.4) 

(3.5) 


<p n (z)<p n+ i{z)w(B) d0 

2 7T _ 

<p n (z)z 2 <p n+1 (z)w(9) d9. 

The first integral on the right is zero because of orthogonality. For the second integral we 
use the recurrence (3.2) for the orthonormal polynomials 




Zipnfz) 


~^(p n+1 (z) +a n ip* n (z ) 

^n+1 


(3.6) 
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to find 


A 

Air 



z 2 (p n (z)(p n+l (z)w(Q ) do 


K n X 


-2?r 


z(p n+1 (z)ip n+1 (z)w(9) dO 


K n+ i Air J 0 

y r 2?r_ 

+ an Air / Z( Pn( z ) ( Pn+i(z)w(9)d9. 


If we use (3.6) for n + 1 and orthogonality, then 


K n X 
K n +l Air 


-2?r 


^ n+ i(^)93 n+ i(^)w(6») d0 


because 


®-n +1 


Kr, 


X 


-2?r 


K n+ i 47 t y 0 


¥>n+l(*Vn+l(*M0) 


A _ 

0 ^n+l^n 

^ ^n+1 


^n+l(-) 


^n+l(0) 

K'n+l 


(/? n+ i(z) + lower degree terms 


—ck„ </? n+ i(z ) + lower degree terms. 


In a similar way we have 

z<p*(z) = —oin-i —— ip n+ \(z) + lower degree terms 

^n+l 


so that 


a r 


X_ 

1 Air 


-2tt 


A 


Kn 


zip^(z)y n+1 (z)w(Q) dd = ——OL n a n -i- 

2 K n +l 


Combining all these results gives 


dz X K r 


-L- I , . / \ S' // \ ^ 1 ^Tl / _ > 

— / (p n (z)(p n+ i{z)w (z) — = -- {a n+ 1 a n + a n a n - 1 ). 


2m 


z 2 K n+ 1 


(3.7) 


We can compute this integral also using integration by parts, to find 


dz 


— / (p n (z)(p n+1 (z)w\z) — = 


2iri 


2iri 


[zip n (z)ip n+l (z)]'w(z) dz. 


We have to be a little bit careful because z<p n {z) is not analytic in the complex plane, but 
on the unit circle T we have <p n (z) = z~ n <p*(z) so that 


1 


— / [z<p n (z)<p n+1 (z)]'w(z) dz = — I [z n 1 (f* n (z)(f n+1 (z)]'w(z) dz 


2iri 
n + 1 


'T 


I* rj y 

/ Z~ n -Vn(zW n+ l(z)w(z) - 
J T z 

+ f - 


2iri 


+ f z n <Pn( z Wn +1 (*)«>(*) y- 
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If we use the recurrence relation (3.6) then 


Z U Vn(*Vn+l(*)«K*) ~ = 2 ^ / -^n(-)^n+l(-)w(0)rf0 


by orthogonality we find 


* n [^n]'(-)^n+l(-)w(^) 


Z 27T 

= 0, 


z(p' n (z)(p n+1 (z)w{9) dd 


and if we use 


tp' n+ Az) = (n + 1 ) — tp n (z) + lower degree terms, 
ftn 


* n <Pn(z)<Pn+l( Z )™( Z ) 


z 2n 


ip n (z)<p' n+1 (z)w(9) d0 


= (n + 1 )- 


These computations give 


— / [zip n (z)y n+1 (z)]w(z) dz = (n + 1 )- 

Z7TI Jj 

Now we can combine (3.7) and (3.8) to find 


A ft n / - . -\ _ / . \ K n +1 (^ K n \ 

0 (^n+l^n T C^nC^n—l) (ft "hi) 1 1 o ) 

2 ftn+1 ftn \ ftn+1 / 

which, together with (3.3) gives 

~(1 |c^n| ) (c^n+l^n "P QC n CV n — l) (ft “1“ l)|c^ n | . 

Recall that w(—6) = w{6) implies that the a n are real. Hence when a n ^ 0 then 

— — (1 — a^)(a n+ 1 + a. n - 1 ) = (n + l)a n . 


This non-linear recurrence relation corresponds to the discrete Painleve equation d-Pn 

x n (an + (3) + 7 
37,+1 T 2-n—1 1 n 

with a n = x n , a — (3 — —2/A and 7 = 0. The initial values are 


— —1, OIq — 


f o 2n zw(0 ) dd = /i(A) 

f^w(9)d6 /o(A) 
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4 Discrete orthogonal polynomials 

In this section we will study certain discrete orthogonal polynomials on the integers N. 
The orthonormality now becomes 

OO 

^Pn(k)p m (k)w k = 5 n , m , n,m> 0. (4.1) 

k =o 

Instead of the differential operator we will now be using difference operators, namely the 
forward difference A and the backward difference V for which 

A/(x) = f(x + 1) - f(x), V/(x) = f(x) - f(x - 1). 

We now have two sequences {a n : n = 1,2,...} and {b n : n = 0,1, 2,...} of recurrence 
coefficients, and we need two recurrence relations to determine all a n and b n . 


4.1 Charlier polynomials 


Charlier polynomials are the orthonormal polynomials for the Poisson distribution 

n k 


Wk = 


kV 


fceN, a > 0. 


Observe that 


w k -i = - w k (4.2) 

a 

which is the (discrete) Pearson equation for the Poisson distribution. It can also be written 
as aVwk = (a — k)w k - The Pearson equation gives the following structure relation for 
Charlier polynomials. 

Lemma 4.1. For the orthonormal Charlier polynomials one has 

Pn(x + 1) = p n (x) + —p n - l(x), (4.3) 

a 

where a n are the coefficients in the recurrence relation (1.2). 

Proof. If we expand p n {x + 1) into a Fourier series, then 


p n (x + 1) = ^ A n jpj(x), 

3=0 

and if we compare the leading coefficients then A n;n = 1. The other Fourier coefficients 
are given by 

OO OO 

A n>j = ^ p n {k + 1 )pj(k)w k = J2Pn(k)Pj(k - l)w k -i, 

k=0 k=1 

and if we use (4.2) then this gives 


A n ,j = ~ ^2 kp n (k)pj(k - 1 )w k . 


k=0 
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The polynomial xpj(x — 1) has degree j + 1, hence by orthogonality A n j = 0 whenever 
j < n — 1. For j — n — 1 we have 

xp n -i(x — 1) = ——— -pn(x) + lower degree terms 

In 

so that (1.4) gives the desired result. □ 

Note that we can write (4.3) also as 

A p n (x) = —p n -i(x). 
a 

If we compare the leading coefficient in the latter, then nj n = a n 7 n _i/a, so that 

= an. 

We will now compute the sum 

OO 

22 p ^ k + l)w k 

k= 0 

in two different ways. If we use the Pearson equation (4.2) then 

OO OO 1 OO , 

22 p ^ k + 1 )w k = ^p 2 n {k)w k _ i = -^2 kpl(k)w k = -A. 

LL CL 

k= 0 k= 1 k= 0 

On the other hand, the structure relation (4.3) gives 

oo oo 2 2 

22Pn( k + l ) w k = 22 ij )n (k) + C ^p n -iik2j w k = 1 -h—. 

k =0 /c=0 

Combining both computations gives 

a n 

b n = o H-- = n + a. 

a 

These simple computations show that the recurrence coefficients for Charlier polynomials 
are given by 

a n = v^an, b n = n + a. 


4.2 Generalized Charlier polynomials 


If we take the weights 


w k 



k G N, a > 0, 


with N G {1, 2, 3,...}, then for IV = 1 one finds the Charlier polynomials and for N > 2 
the generalized Charlier polynomials. These were introduced by Hounkonnou et al. in 
[10]. The Pearson equation is 

k N 

w k -i = — w k , (4.4) 

a 

which can also be written as aS7w k = (a — k N )w k . For N > 2 the factor a — k N is a 
polynomial in k of degree greater than one, and hence the weight is no longer classical 
but semi-classical. We will investigate the case IV = 2 in more detail. 
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Lemma 4.2. For N = 2 the generalized Charlier polynomials satisfy the structure relation 

Pn(x + 1) = p n (x) + —Pn-l(x) + ^-Fzl p n _ 2 (x), (4.5) 

a n a 

where a n are the recurrence coefficients in the three-term recurrence relation (1.2). 

Proof. If we expand p n {x + 1) into a Fourier series, then 

n 

Pn(x+ 1) = A nJ Pj(x). 

j= 0 

Comparing coefficients of x n gives A n;n = 1, and comparing coefficients of rc n_1 gives 
A n . :n -i = n/a n . The remaining Fourier coefficients are given by 

OO OO 

And = T.Pnik + 1 )pj(k)w k = J2p n (k)p 3 (k ~ l)^fe-l- 
k= 0 fc=l 

If we use the Pearson equation (4.4) then 

^ OO 

A n j = - y^P n (k)pj(k - 1 )k 2 w k . 

a z —' 

k=0 

The polynomial x 2 pj(x — 1) is of degree j + 2 and hence by orthogonality = 0 for 
j < n — 2. For j — n — 2 we have 

x 2 p n - 2 (x — 1) = ~ p n {x) + lower degree terms 


so that 


A n .n —2 — 


1 'Tn—2 Q'nQ'n—1 


where we used (1.4), which gives the required result. 
The structure relation (4.5) can also be written as 


A Pn (x) = - Pn ^(x) + ?A^zl pn _ 2 ( x) ' 
a„ a 


If we compare coefficients of x n 2 , where we use 

p n (x) = 7 n x n + 8 n x n ~ l + 


then we hnd 


/ . N, r. 77/ _ C^T) C^T) -1 

7 n + ( n ~~ l)^n = - 8 n -l H-7n—2- 

CLft Oj 


If Xi tn < X 2 ,n < • • • < x n ^ n are the zeros of p n , then by Viete’s symmetric formulas we 
have 

r- n 

On \ A 

/ %k,n • 

'v ^^ 
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The zeros of p n are equal to the eigenvalues of the truncated Jacobi matrix 


/&0 





w 

Cl2 



Cl2 

b 2 

03 



a 3 

b 3 


\ 


V 



1 I 

bji— 1 / 


and the sum of all eigenvalues is the trace of the matrix, hence 


In 


n— 1 

Y. h r 

3=0 


If we use this in (4.6), then 


nb n -i + b 

3=0 


n —1 o 9 

_ «-1 
'3 


a 


(4.7) 


In order to get rid of the non-homogeneous terms, we put b n — n + b n , and the relation 
becomes 

n —1 9 9 

i , ffiX-l 

-n6 n _i + bj =-. 

i=o a 

Differencing both sides gives 

-na(b n - b n - 1 ) = a 2 (a 2 +1 - a 2 ^), (4.8) 

which may be considered as the first Freud equation for the recurrence coefficients. 

Next, we will compute 

OO 

y^/Pnjk + l)p n -i(k + 1 )w k 

k =0 

in two different ways. First we use the Pearson equation (4.4) to find 


^p n {k + l)p n -i(k + 1 )w k 

k =0 


y^Pn(fc)Pn-l(fc)Wfc-l 


k =1 


y^/c 2 p w (A:)p w _i(fc)w fc 


fc =0 


-(J 2 ) nl „_l. 


The entry (J 2 )n in _i can be computed easily by repeatedly using the recurrence relation 
(1.2) and is equal to (J) 2 n _ x = a n {b n + b n - 1 ), so that 


OO 

^2p n (k + l)p n -i(k + l)w k 

k =0 


Q n(p n ~b b n — i) 
a 


(4.9) 
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On the other hand, we can use the structure relation (4.5) to find 

OO 

y ^Pn(k + l)p n -i(k + 1 )w k = 

k =0 


where the last equality follows from the orthonormality (4.1). Combining (4.9) and (4.10), 
and recalling that b n — n + b n , then gives the second Freud equation 

a 2 n (b n + b n - 1 + n ) = na. (4.11) 

If we eliminate na from the two equation (4.8) and (4.11), then 

~ib n + 6 n _i + n)(6 n - b n -\) = a 2 n+l - a 2 n _ v 

Summing both sides of this equation gives 

n— 1 

b 2 n + Y bk ~ nb n = a 2 +1 + a 2 n - a. (4.12) 

k =0 

Summing both sides of (4.8) gives 

™-i \ 

E h k ~ U bn I = a n a n+l ■ ( 4 - 13 ) 

k =0 / 

Combining (4.12) and (4.13) then gives 

a n a n+ 1 — a ( a n +1 + °n) + ^ = 


which is equivalent to 


ab 2 n = (a 2 +1 - a)(a 2 l - a). 

(4.14) 

This means that — a and a ^ +1 — a have the same sign, and since a 0 

conclude that a 2 — a < 0 for n > 0. We may therefore write 

= 0 we must 

o(l Ci), 

(4.15) 

with Co = 1, and then (4.14) becomes 


bn ^ n (‘n ■ 1 ■ 

(4.16) 

The second Freud equation (4.11) becomes 


(1 Ci)v^®(Ci+i T Cn—i) nc n . 

(4.17) 


If we compute the coefficients c n from the recurrence relation (4.17), then we obtain the 
recurrence coefficients b n — n + b n from (4.16) and the a n from (4.15). The non-linear 
recurrence relation (4.17) corresponds to the discrete Painleve 11 equation 

x n {an + (3) + 7 
Xn +1 T 1 9 

1 - 4 



OO / 

E \Pn( k ) + —Pn-l(k) + 2(fe) 

V 

x (p n _i(A;) + -—-p n - 2 (&) + Pn _ 3 (k )^ 

^ + (n ~ 1)a ", (4.10) 

a n a 
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with c n = x n and a = 1 /y/a, (3 = 7 = 0. We need to find the solution with c 0 = 1 and 
c\ = 1 — a 2 /a. Observe that if we require that pi(x) — (x — &o)Po/°i is orthogonal to p 0 , 
then 

OO 

- b 0 )w k = 0 

k =0 


so that 


where 


b 0 = 


E k 

k =0 


w 


00 


£ 



y/ah(2y/a) 

/ 0 (2\/a) 


z ) = 

k =0 


(z/2) 2k+u 

k\Y(k + v + 1) 


is the modified Bessel function. From (4.16) we then see that 


/i(2\/a) 

J 0 (2\/a) 


(4.18) 


The non-linear recurrence relation (4.17) with initial conditions Cq = 1 and Ci = 
/i( 2 y / a)/J 0 ( 2 v /a) is again very unstable for computing all the c n recursively. One can 
show that the discrete Painleve equation with 7 = 0 and Co = ±1 has only one solution 
for which — 1 < c n < 1 for all n > 1 (see [24]), and this is the solution that we need 
since a 2 n = a(l — c^) needs to be positive. Hence a slight deviation from the actual initial 
value Ci will destroy the positivity of the o? n eventually. In Figure 2 we have plotted the 
c n obtained from the recurrence relation with an accuracy of 30 digits. The c n converge 
quickly to zero, but for n near 40 we see that the c n deviate quite a lot from zero. 

The discrete Painleve 11 equation also arose in Section 3 for the Verblunsky coefficients 
of certain orthogonal polynomials on the unit circle. Verblunsky coefficients always have 
the property that |a n | < 1 for n > 0 , hence in that case one also requires the unique 
solution with a_i = —1 for which — 1 < a n < 1 for n > 0. Observe that there is a shift in 
the index since we are using Verblunsky coefficients, in which case the recurrence starts 
with a_i = — 1 . 

Obviously the equation (4.17) satisfies the discrete Painleve property. Indeed, we have 


Cn+1 


nc r 


Cn—!•> 


hence a singularity will appear in c n+ 1 whenever c n = ±1. A careful analysis gives that 
for c n near 1 


Cn 

C'n ■ 1 


Cn+2 
Cn -(-3 


1 + e 


n 


2 y/a e 
-i + n + 2 


n 


77 

—= - c n _ 1 + O(e) 
4 y/a 

e + 0(e 2 ) 


71+1 

+ 2 ) 


n 

77 + 2 


Cn-1 + C*( e )i 
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30 digits accuracy 


: 



O 

: 



« 

_ 0 0 Q 

0 0 0 0 0 0 



O 

10 

0 °20' ° 

' ° °30° ° 

°40° 6 0 50 




O 




° 




O 


Figure 2: The result of computing c n from d-P n using 30 significant digits 


and near — 1 


c n — — 1 + e 

77 1 77 

cn + i = - TT - + TT -c n - 1 + 0(e) 

2 y/a e 4 y/a 

c n + 2 = 1 + e + 0(e 2 ) 
n 

n + 1 n , . 

c n +3 — -7=7 —-777 -—7 c n-i + 0(e), 

yja{n + 2) n + 2 

so that in both cases the singularity is confined to c n+ 1 and c n+ 2 . Observe that the critical 
value 1 for c n results in the critical value —1 for c n+ 2 , and that the critical value —1 for 
c n results in the critical value 1 for c n+2 . 


5 q -Orthogonal polynomials 

Here we consider orthogonal polynomials on the exponential lattice {±q n , n G N}, where 
0 < q < 1. The orthogonality is of the form 


Pn(x)Pm(x)w(x) d q X = 6 m , n , 


(5.1) 


'-1 


where the (/-integral is defined by 

/: 


f(x) d,x = (l-g)J2 f(l k )q k + (1 - q) J2 


k =0 


k =0 
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We will only consider even weights for which w(—x) = w(x), in which case the orthogonal 
polynomials have the symmetry property p n {—x) = (—1 ) n p n (x), i.e., the polynomials are 
even when n is even and odd when n is odd. The recurrence relation will then be of the 
form 

xp n (x) = a n+l p n+l (x) + a n p n - i(x), (5.2) 

with i = 0. The results in this section were obtained for the first time by Nijhoff [19], 
but we take a slightly different approach. 


5.1 Discrete (/-Hermite I polynomials 

The orthonormal discrete (/-Hermite I polynomials [12, §3.28] are given by 


Pn(x)p m (x)(qx ; q)oo(-qx ; q )oo d q x = 8, 


771 , 71 -) 


1 


where 

OO 

(x;q) oo = JJ(1 ~xq k ). 

k =0 

Observe that 

OO 

(qz; q)oo(-qx ; q) oo = JJ(1 - x 2 q 2k ) = ( x 2 q 2 ; <? 2 )oo, 

k =0 

so that the weight can be defined as w(x) = (x 2 q 2 ] q 2 )^- In terms of the ^-exponential 
function E q (z) = (—z;q) 0 0 we have w(x) = E q 2 (—x 2 q 2 ), and since E q (( 1 — q)z) — * exp(^) 
when q —> 1 it follows that w(^/l — q 2 x) —> exp(— x 2 ) when q —> 1, which shows that this 
weight is a (/-analog of the Hermite weight. One easily finds that 


(1 — x 2 )w(x) = w(x/q), 


(5.3) 


which is the Pearson equation for this weight on the (/-lattice. The structure relation for 
the corresponding orthogonal polynomials is in terms of the (/-difference operator D q for 
which 


D q f(x) 


f(qx) - f(x) 
x(q — 1) 

m, 


if x 7 ^ 0, 

if x = 0. 


Lemma 5.1. The discrete q-Hermite I polynomials satisfy 

D q p n (x) = ^ ‘ Pn ~ 1 - 


(5.4) 


Proof. Clearly D q p n (x) is a polynomial of degree n — 1 and D q p n (—x ) = (—l) n 1 D q p n (x). 
If we expand this polynomial into a Fourier series, then 


71—1 

D q p n (x) = ^ apnPjjx ), 
1=0 


with 


-l 


dj,n = / D q p n (x) Pj(x)w(x) d q x. 
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The symmetry shows that a j n = 0 whenever n — j is even. When n — j is odd then 

OO 

aj, n = 2 (1 ~ q)^2D q p n (q k ) Pj (q k )q k w(q k ) 

k=0 
oo 

- -2 - Pn(l i ))pA<lX‘l k ) 

k=0 

oo oo 

= -2 +2^R,(9*)pj(9*)to(9‘). 

k=0 fc=0 


Both sums are finite since either p n or pj is an odd polynomial. Using the Pearson equation 
(5.3), and a shift in the summation index in the first sum, gives 




-2 - Q 2k ) + 2Y,Pn(q k )Pi(q k )w(q k 


k=0 


k=0 


Pj{x/q) -Pj(x 


p n (x r n "c w (x) d q x + - 

x(q — 1 ) 1 


xp n (x)pj(x/q)w(x) dqX. 


mi x(q- L) ' ■ i -qj -1 

The first integral on the right is zero because of orthogonality. The second integral only 
gives a contribution when j — n — 1, in which case 


^71—1,71 


xp n {x)p n -\{x/q)w(x) d q x. 


1-9 J -1 

The recurrence relation (5.2) gives 

xpn-xix/q) = qa n p n (x/q ) + qa n -ip n - 2 (x / q ), 

and since 

Pn{x/q ) = q~ n p n {x ) + lower degree terms 

this gives 


®n—l,n 




1 - q J- 

which gives the desired structure relation 


Pn(x)p n (x/q)w(x) d q x = ™ 

l 9 pl — o) 


□ 


If we compare the leading coefficients on both sides of (5.4), then 

1 - q n _ a n 

7n 1 -q ~ 7n_1 q n ~ l [l — q) ’ 

so that we find 

a 2 n = q n -\l-q n ), 

which are indeed the recurrence coefficients as given in [12, §3.28]. So for these orthogo¬ 
nal polynomials the recurrence coefficients can be found immediately from the structure 
relation (5.4). Observe that the o 2 n tend to zero exponentially fast and that 

lim a 2 Jq n ~ l = 1, 


and 


lim 

i 1 — q 2 


n 

2’ 


and the latter are the recurrence coefficients (2.8) for the Hcrmite polynomials (p = 0). 
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5.2 Discrete g-Freud polynomials 

A (/-analog of the Freud polynomials with weight exp(— x 4 ) can be obtained by taking 
the weight w(x) = (g 4 x 4 ; g 4 ) 0O = E q i{—q A x A ) on the exponential lattice. Observe that 
w (W q 4 x) —> exp(— x 4 ) as q —> 1. This weight satisfies 

(1 — x A )w(x) = w(x/q), (5.5) 


and the structure relation for these semi-classical polynomials is: 

Lemma 5.2. The orthonormal polynomials for which 

r± 


satisfy 


with 


Pn(x)Pm(x)(x q \q ) OQ d q x = < 5 m , n , 


-1 


D q p n {x ) = 


— 


p n -l(x) 

+ 1 

1 ~q 

i -q 

^n^n—l^n—2 


q n-3 


(n +1 

n —2 

a n ( ^2 

q n ~i l / 

y \j=1 


3 = 1 


Pn-3(x), 


■ 


(5.6) 


(5.7) 

(5.8) 


Proof. Expanding D q p n into a Fourier series gives 


n—1 

DqPn(x) = YapnPAx), 

j =0 


with 


-i 


Uj,n = / D q p n (x) Pj(x)w(x) d q X. 


Again « jn = 0 whenever n — j is even. When n — j is odd then, as in the proof of 
Lemma 5.1, we have 


® j,n 


Pn\X 


'-1 


^ Pjfx/q)-Pj(x) 
x(q — 1) 


w(x) d q x + 


1 -q J- 


x i p n (x)pj(x/q)w{x) d q x, 


where we have now used the Pearson equation (5.5). Again the first integral on the right 
vanishes because of orthogonality. The second integral only gives a contribution when 
j — n — 1 or j — n — 3. For j — n — 3 we have 


3,n 


QJ- i 


x 3 p n (x)p n _ 3 (x)w(x) d q x, 


and since 


x 3 p n - 3 (x/q) 


|/ " 3 q n+3 p n (x) + lower degree terms 
Kn 
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we easily find 


tin— 3,n 

which gives (5.7). For j = n — 1 we have 

1 f 1 


^n^n—l^n—2 

(1 — g)g n_3 ’ 


^71 — 1,71 


x 3 p n (x)p n _ l (x)w(x ) 


and if we write 

x 3 p n _ 1 (x/q ) = A n+2 p n+2 (x) + B n p n (x) + lower degree terms, 
then the orthonormality gives 


B n 


®n—l,n 


1 - g 


If we compare coefficients of x n in (5.9) then 

fin-iq " +J> — A n + 2 $n +2 + B n 7 n , 
where p n (x) = 7 n x n + S n x n ~ 2 + • • •, so that using (5.7) gives 


t-, Q'n ( 2 fin— 1 ^n+2 
Br, = -- I q - 


<?” 1 V" 7n-l 7n+2 
If we compare coefficients of x n_1 in the recurrence relation (5.2) then 

fin ®n+lOn+l “h &n~ln— 1; 

from which one easily finds 

On+l 0 n _ ^ 2 


7n+l 7n 


which gives 


fin 

7n 


71—1 


£7 


i=i 


and using this in the formula for B n gives the desired expression (5.8). 

If we compare coefficients of x n in the structure relation (5.6) then 

7n(l - Q n ) = ln-\B n . 

Comparing coefficients of x n ~ 3 in (5.6) gives 

^n(l Q U 2 ) = Sn-\B n + 7 n sA n , 

which together with (5.11) gives 


An - 


7n (fin,* n -2\ fin-1 


(1 _ q n->) _ _2Zi(l - g* 

'"'In —3 VTti 'TVi—1 


(5.9) 


(5.10) 

□ 

(5.11) 
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Together with (5.7) and (5.10) this gives 


aX-^n-2 = Q n 3 1 2 


n—2 

9" i 1 - <? 2 ) X a i “ ( X “ ^' 2 ) a l-i 

3= 1 


(5.12) 


On the other hand, if we compare (5.11) with (5.8) then we find 

/ n -\-1 n—2 

r „ . 

3 


1 - q r 


a r 


CL r . 


T 


T 5>,-<? 2 Va 2 


d =i 


X u i / ’ 

1=1 


which can be written as 


n—2 


q n 1 (! - 9 n ) = a n+i + « 2 + a 2 -i + (! - 9 2 ) X 1 


(5.13) 


l=i 


If we take (5.12) with the index n raised by one, then we can find 


n—2 


(! - q 2 ) Xs 2 = ? 2 n+ 3 0 2 +l 0 2 «n-l - (! - 9 2 )«n-l - (1 - q n+ Vn 1 
1=1 

and if we insert this in (5.13) then we find the second order non-linear equation 

<f~\ 1 - ?”) = <4 (<* 2 +i + + 9 2 aLi + 9 _2 ” +3 “i+i“i;“5-i) ■ 


(5.14) 


We claim that this equation is a (/-deformation of the discrete Painleve I equation. Indeed, 
if we take x n = a 2 / yjl — q 4 then 

1 — q n 

q n ~ l - _ 4 = (x„+i + <T n+1 z n + q 2 x n _i + (1 - g 4 )g“ 2n+3 x n+ ix n x n _i) , 

which for q — > 1 converges to 

U , , , V 

^ *^n(*^'n+l T X n T 2) n _i), 

which is the discrete Painleve I equation (2.12) for Freud polynomials (with p = 0). If we 
put a 2 = q n ~ 1 y n then (5.14) can be rewritten as 


i n (Vn+iyn + l)(</»-li/n + 1) = 1 - 111- 


(5.15) 


This could therefore be called a (/-discrete Painleve I equation (q-Pi). 

We can easily find the asymptotic behavior as n —> oo. First observe that from (5.14) 
we find the upper bound 

<f" + X<<r‘(i-4 n ). 

so that a 4 < q 2n ~ 2 (l — q n ), and a n tends to zero as n —> oo. Let A = limsup^^ a 2 /g n_1 , 
then if we take a such that a 2 /?” -1 converges to A, equation (5.14) gives A 2 = 1. A 
similar reasoning also shows that B = lim inf„^oo a^ l /q n ~ 1 is such that B 2 = 1. Hence we 
may conclude that 

lim a 2 Jq n - 1 = 1. 

n—xx) 
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The equation (5.15) has the singularity confinement property. Indeed, a singularity 
occurs for y n+ \ whenever y n = 0. So if we put y n = e, then some straightforward calculus 
gives 


Un 
Un +1 

Dn +2 

Dn +3 

Un +4 


e, 


<?-"(! - 9”) - - <TV-l + C( £ ), 

“ + Vn-i/q + e(e), 

-V 2i r^— e + 0(6 2 ), 
l-g n v ' 

2 i-g w 

1 - q n + 3 


q 


Un -1 + C*( e )- 


Hence the singularity is confined to ?/ n +i, y n + 2 , Un+ 3 - 

Again the recurrence relation (5.14) or (5.15) is very unstable for computing the re¬ 
currence coefficients recursively. One can show [24] that there is again a unique solution 
of (5.15) with yo = 0 which is positive for all n > 0, and this is the solution for which 
y n = a^ l /q n ~ 1 . This solution is such that y n —> 1 and 


V\ = a\ = 


fix x 2 (x 4 q 4 ; q 4 )oo d q x 
J^(x 4 q 4 ; g 4 )oo d q x 


(g; g 4 )oo 
(g 3 ;g 4 )oc 


where the integrals can be computed using the (/-binomial theorem. In Figure 3 we have 
computed log \y n \ recursively for q = 0.9 with 50 significant digits. 



50 digits accuracy 

4- 

O 


O 

3- 

• 

2- 

o 

1- 

• 


0 

0. 

-1- 

10 o 0 0 „ o20o O o o o o ° O o^o » ° o ° “ “ “ u 40“ “ 50 

o ° 

O 

O 


Figure 3: The result of computing log \y n \ from (5.15) (q = 0.9) using 50 significant digits 
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5.3 Another discrete (/-Freud case 

If we take the weight w(x) = (x 2 q 2 ] q 2 )oo(cx 2 q 2 ] q 2 )oo, with c < 1, then w is positive on 
the (/-lattice and it satisfies the Pearson equation 


(1 — x 2 )(l — cx 2 )w(x) = w(x/q). (5.16) 

If c = —1 + a a/1 — q 4 then w(f/l — q 4 x) —> exp(— x 4 — 2 ax 2 ) so that this gives us a q- 
deformation of the Freud weight exp (—a; 4 — 2 ax 2 ). Observe that c = — 1 gives the discrete 
(/-Freud polynomials considered in the previous section and c = 0 gives the discrete q- 
Hermite I polynomials. 


Lemma 5.3. The structure relation for the orthonormal q-polynomials with weight w(x) = 
(x 2 q 2 ]q 2 ) 00 (cx 2 q 2 -,q 2 ) 00 on the q-lattice {±q n , n G N} is 


where 


D q Pn(x) = -^—p n - i(x) + —^ —p n _ 3 (x), 

1 — q 1 — q 


A 

B 


^n^n—l^n—2 



n —2 




(5.17) 


(5.18) 

(5.19) 


Proof. The proof is a straightforward copy of the proof of Lemma 5.2, except that one 
uses the Pearson equation (5.16). The Pearson equation contains the quartic polynomial 
(1 — x 2 )(l — cx 2 ) = 1 — (1 + c)x 2 + x 4 so that one ends up with integrals of the form 


A r , 


B n 



CX 3 ]p n (x)p n - 3 (x)w(x) d q X 
cx 3 ]p n (x)p n _ 1 (x)w(x) d q x 


cA n 


—cB n + (1 + c)- 


1 n —1 : 


where A n and B n are given by (5.7)-(5.8). □ 

Reasoning in the same way as in the previous section, i.e., comparing the coefficients 
of x n_1 and x n ~ 3 in (5.6), one arrives at 


q n ~\l-q n ) = 


~ Ca l ( «n+l + Q n+1 «n + ^<-1 ~ 


J2 „2 


1 + C 


— 2n+3 2 22 

-eq + 


(5.20) 


If we put a 2 = q n l y n , then this can be rewritten as 

(1 - y n ){ 1 - cy n ) = q n (cy n+1 y n - l)(cj/ n _ij/ n - 1), (5.21) 

which is a more general form of the (/-discrete Painleve I equation in (5.15). 
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Appendix 

Several discrete Painleve equations have appeared in the literature, and the list is certainly 
longer than the six Painleve differential equations. Sakai [ 21 ] made a classification in terms 
of rational surfaces associated with affine root systems and the most general (elliptic) 
discrete Painleve equation is related with the affine Weyl group symmetry of type E 8 . 
Sakai’s classification does not give explicit expressions for the discrete Painleve equations. 
A few important discrete Painleve equations are listed below. The list was compiled by 
Peter Clarkson and I thank him for his permission to present it in this paper. 


A.l Discrete Painleve equations 

1 O , Z n + 7(-l) n , x 

d-Pi x n+ i +x n + x n _i = -b o 


x r 


d-P n x n+l + x n _i = 


x n z n + 7 

1-4 


d-Piv (x n+1 + x n )(x n + x n _i) = 

d-P v 


(4 - 4) (4 -/i 2 ) 


(x n + z n y - y 2 

4n+l T X n ^n+1 %n) (x n T X n —\ Z n Z n — \) [(^n Z n )~ ® ] [4n %n ) ft ] 


(*^n +1 T X n )(x n T X n —\) 
where z n = an + ft and a, ft, 7 , 5, k, fi are constants. 


(x n - 7 2 )(x n - 5 2 ) 


A.2 (/-discrete Painleve equations 


<?-Pll (■ X n +lX n - l)(x n x n _i - 1) = 


^nftn—lXn 
Qt (x n aX n ) 


Q Pi I 2'n+l‘^'n—1 ®A 

(/“Pill X n -\-\X n —1 


X n T X n 

x n (x n 1 ) 

(x n + a){x n + ft) 


q- Piv (■ x n+1 x n - i)(x„x n _i - 1 ) = 


q- Py {Xn+lXn - 1) {x n X n —l ~ 1) = 


ft)XnX n T l)(dX n x n A 1) 

7 ${x n + o/){x n + 1 / a)(x n + ft)(x n + 1/ft) 


( / ~)X n x n T l)(5A n x n T 1) 

'fSXKxn - a)(x n - 1 / a)(x n - ft)(x n - 1/ft) 


( x n r yX n )(^x n SX n ) 


q -P 


VI 


(^n^n+l A n A n _)_i) (x n X n —i X n X n -i) 


(x n x n +1 l)(x n X n —\ 1) 

_ (.Xn (%X n )(x n X n /a) {x n ftX n )(x n X n /ft) 

(x n - y)(x n - 1 / 7 ) (x n - S)(x n - 1/5) 
with a + b + c + d = 0 , p + q + r + s — 0 


where A„ = Ao</ n and a, ft, 7 and A are constants. 



A.3 Asymmetric discrete Painleve equations 


tt-d-Pi 

tt-d-Pn 

a-d-Pm 

a-d-Piv 

a-d-Py 


a-q- Py 

a-d-Pyi 


fl(y n ) + x n f 2 (y n ) 

Xn+1 h{yn) + x n f 4 (y n ) 


9l(x n ) +Vn9 2(x n ) 
9s{^n) "P Und^illn) 


*£n+1 P %n P /Jn d P 


^n ~ 7 

Un 


y n + y n _i + = 5 + 


x ri 


*^n+1 “P *^n 


Un P Vn—l — 


2 (jjn^n P T) 

i-yS 

2(^n^n+l/2 

1 — xi 


%n+l%n 


(y n - g n a)(yn - g n &) 


UnUn —1 


(x n — q n a)(x n — q n /3) fa/3 ab 

^ = q ffd 


(Vn P *^n) (*^n+l P Vn) 
{/Jn P 2'n)(*^n P Vn—l) 


(y n ~c)(y n -d) (x n - 7)(x n - 5) 

(y n - a) (y« - b ) (yn - c) (y„ - d) 


(yn + 7“'Zn)(yn-7—'2 ; n) 

(x n + a)(x n + 6)(x n + c)(x n + d) 


(x n P d Zn+l/ 2 ') d %n+ 1/2) 

with a+6+c+d=0 


(JJn P *^n+l ^n+l) (*^n+l P Vn Z n Z n — l) 

{.Vn P 3p)(*^n+l P Vn) 

_ (y n - z n - a) (y n - z n - &)(?/„, - z n - c)(yn - z n -d) 

(■Vn ~ P) (Vn - 9) (yn “ 0 (?/„ - s) 

(i yn P *^n+l ^n—1/2) (*^n+l P Vn ^n—1/2 ^n—l) 

{/jn P Xji){x n P yn— l) 

(*^n ^n—1/2 P o)(Xjj Z n —\/2 P b)(x n Zn—1/2 P C)(x n ^n—1/2 P d) 

(^n p p) (P n p g) (x n P r) (x n + s) 


(x n i/ n - i)(x n _i2/„ -1) = q 


2 n 


(yn - a)(y n - 6)(y n - c)(y n - d) 

(g n - ny n )(q n - y n /n) 


(x n y n - l){x n y n+ i - 1) = g 2n+1 — 
with abed = 0 


l/a)(^n - l/fr)QEn - l/c)(^n ~ 1/d) 
( g n+l/2 _ ^y n )(q n + 1/2 _ ^//j) 


{ ^n^n+l 
HuUn-l 


fePijyn - g^lKPn - g n «2) 

(//n ){/Jn ^ 4 ) 

«3«4(^n ~ ~ g^Ai) 

(^n - /5 3 )(x„ - P 4 ) 


with 


a±a2 _ P1P2 

Oi%a 4 P 3 P 4 


29 



where z n = an + (3 and a, b, c, d, p, q, r, s, a, (3, 7 and <5 are constants. 


A.4 Alternative discrete Painleve equations 

,-n , z n + 7( — 1)” , 

a-d-Pi x n+ i + x n + x n _i =-h /i 


1 2 
H-= —x; 


*^n+l ~t~ X n ~h X n _i 

, _ Z n , 7 

372+1 + Xn-1 — -1- 

x n x n 


Z n 

x n+l ~P X n — i 

•En 


x n+l x n —1 


exp(^n) 7 

' 9 

IJQ 


a d Pjj x n ^\ H“ 


+ 


3^ n + 7 


x; 


2ra-l 


a^n+ia^n "P 1 x n x n —\ ~{~ 1 


~ x n 3 “t" ^ri 


a-d-P 


v 


(a-n+i 3 ~ x n ‘ 2 z ri ^ i ^(^ x n *h x n —\ (372 ^22+1/2) ^ 


(a-n+i “h x n )(x n -P x n _i) 


^n+ 1/2 


+ 


Zn—l/2 


1 X n X n -|_i 1 X n X n _i 


— I 1 + Z n + 


(®n - 7) 2 

KX„ 1 — X 


+ 


22 M 


(1 + X n ) 2 1 + X, 


[I* + 


where = an + (3 and a, /3, 7 and 5 are constants. 


(-l) n 7] 
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A.5 Other discrete Painleve equations 


i -q . Un^n 4" 7 

d-Pi x n+ i +x n = - 5 -, 

Vn 

DnZn + 7 

vi - A 


d-Pn x n+ \ +x n = 


d-Piv x n x n -1 = 


(y n + z n - a)(y n + z n -b) 


yl~i 2 


with a + b + c + d = 0 


Vn 4" Un—1 

Un 4" Vn—1 
Un + 2 /n+l 


2'n-^n+l/2 4" J 


Xt 


XnZn+ 1/2 4" J 

r 2 — i,2 

Z n + 1/2 + C Ui+1/2 + d 
Xnl +1 X „/7 4 1 


d-Piv x n x n —\ = 


a{y n + Zn - ft) 

y 2 n - 7 2 


i tj , z n -\- fi z n fi 

d-Py a; n + x n _i = —-- + , 

1 + y n /t 1 + ty n 

. _ 7 , + <5 

4" 1 4" > 

1 “1“ Vn 1 2/n 


c Z n+ i /2 + d 
Un + 2/n+l ~-1 


VnUnf-l 




(a: n - . 1 / 2 )' - ft 2 

^2 _ ~,2 

J 'n / 


•^n Az 
VriVn+l 7 9 2 

4 - h 


where z n = an 4 - (3 and a, b, c, d, p, q, r, s, a, (3, 7 , 5, k and p are constants. 
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